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PACS 42 . 50 . Dv - Quantum state engineering and measurements 
PACS 03 . 67 . -a - Quantum information 
PACS 85 . 25 . Cp - Josephson devices 

Abstract. - We introduce a type of measurements that generalize the so-called "partial measure- 
ments" performed in recent years with phase qubits. While in the case of partial measurements it 
has been demonstrated that one could undo the effect of the measurement only for non-switching 
events, we show here that generalized partial measurements can be reversed probabilistically for 
both switching and non-switching events. We calculate the associated Fisher information and 
discuss the estimation sensitivity for quantum tomography. Two ways of implementing this type 
of measurements with superconducting qubits are proposed. 



Introduction. — We consider a type of measurements 
which are generalizations of the "partial measurements" 
demonstrated experimentally in recent years in the field 
of superconducting qubits 1 . An interesting feature of 
partial measurements is that they can be reversed [2] if 
the SQUID used for measurement did not switch in the 
running-wave state. In case the SQUID has switched, 
the qubit is destroyed, and a voltage is registered by the 
measuring device [3], precisely like in quantum optics an 
event would correspond to a photon being absorbed in 
a detector. In contradistinction, our generalized partial 
measurements have the property that they can be proba- 
bilistically reversed for both measurement results (switch 
or non-switch). We further show that such measurements 
can be performed along any direction in space and we cal- 
culate the corresponding Fisher information metric. Fi- 
nally, we give two explicit physical constructions of these 
measurements, one using two qubits and the other using 
a single qubit in a two-well potential. 

Generalized partial measurements: definition. 

— We consider a qubit with states |0) and |1), pre- 
pared in an unknown pure state \ip) = cos(0/2)|O) + 
exp(iip) sin(0/2)|l). Partial measurements Q] are charac- 
terized by a single parameter p, which is the probability of 
switching (or tunneling) from the state |1); tunneling from 
the state |0) is forbidden, and as a result one of the Kraus 
operators defining these measurements is in fact a pro- 
jection. The immediate generalization of partial measure- 
ments is to allow also |0) to switch (tunnel) with probabil- 
ity q. This results in POVM measurements [4] described 



by two measurement operators, M m and Mfn, 



M m = 

Mm = 



^/r^|o)<o| + vT=p|i><i|, a) 

V9|0)(0| + VP|1>(1|, (2) 



and corresponding to two effects E m — M^Mm and 
Em = M^Mm, which are positive operators realizing a 
semispectral resolution of the identity, E m +E r7l = 1. Here 
m and fh ("not m") are measurement results correspond- 
ing respectively to the absence or existence of a switching 
(tunneling) event. As mentioned before, the parameters 
< p < 1 and < q < 1 have the meaning of probabilities 
of obtaining the result fh depending on whether the qubit 
is in the state |1), or |0) respectively: indeed, the probabil- 
ities of obtaining the results m and fh can be immediately 
calculated, 



P„ 
P, 



M-EgV) - l-qcos 2 (#/2)-psin 2 (0/2),(3) 
(V-|£rn|^> = g cos 2 (0/2) +p sin 2 (0/2), (4) 



and the wavefunctions after the measurement are \ip m ) 
and \ipfn), 



--M m \ip) 



(5) 



cos -y/T=q\0) + sin -y/T^pe^\l) 



--M m \ij) 



(6) 



cos -V9|0) + y/pe* sin-|l) 
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Reversal of generalized partial measurements. — 

We now turn to the problem of time-reversing the mea- 
surements: we claim that, if p, q £ {0, 1}, both measure- 
ments defined by Eqs. (JTJ [5]) can be probabilistically re- 
versed, thus bringing back the system to the exact initial 
state. Indeed, if p, q ^ {0,1}, M m and admit an 
inverse, 



m: 



M: 



1 



^{l-p)(l-q) 
1 -XM^X, 



XM m X, (7) 



pq 



(8) 



where here and in the following we will use the standard 
quantum-information notations [5] for the Pauli matrices, 



X 



1 

1 



,Y 



-i 

1 



,Z 



1 

-1 



(9) 

Suppose now we were measuring an arbitrary state 
and the result m has occurred: then, we apply X and 
measure again, and if the result m occurs again, then 
we apply X again and recover \tp)\ if, instead, we get 
the result fh then we fail. The other possibility is that, 
when we first measure \ip), the result fh occurrs: we ap- 
ply X, measure again, and again either fh occurs, in 
which case we apply X and recover \ip)\ or, m occurs, 
in which case we fail. The probability of success in each 
case can be found by using the usual rules of multiplica- 
tion for conditional probabilities at each step or directly 
by (ip\M- 1 XM m X\ip) = (l-p)(l-g) for the path m -> m 
and (iplM^ 1 X MfhX\?p) = pq for the path rh — > fh, giving 
a total success rate of 1 — p — q + 2pq, independent of the 
state. 

One might now wonder if it is possible to re- 
verse the result of the measurement determinis- 
tically: that is, is it possible to reverse also 
the "fail" r esults M ^X M m \\[>) = M m XM m \^) = 
sin(6>/2)e l¥ Vg(l -P)|0) + cos(0/2) y/p{l - q)\l) in the 
scheme above? Up to a global phase factor, any 2x2 uni- 
tary can be written as R z (a)R y (/3)R z ('y), where R y , R z 
are rotations around the axes y, z, R z (a) = exp(—iaZ/2), 
R y (f3) = exp(— i(3Y/2). The reversibility condition is 
R z (a)Ry((3)R z (^)MfnX M m \ip) = after some algebra, 
we find that this implies a — /3 = 7 = (that is, an X 
gate again), and also p = q. This means that the only sit- 
uation when we can have deterministic reversal (reversal 
of all paths) is the trivial case in which the measurement 
operators Eqs. (jll2[) become identity. In this situation 
we also notice that Eqs. (13141) yield P m = P™ = 1/2 
with no state dependence, meaning that no information 
about 9 can be obtained! Of course, M m and being 
the identity effectively signifies that no measurement has 
been done. This proves that it is not possible to reverse 
deterministically a generalized partial measurement. 

Quantum tomography using generalized partial 
measurements: Fisher information. Next, we 



show that generalized partial measurements can be per- 
formed along any direction, once they are experimen- 
tally available along the z direction, as in Eqs. ([IJ [2]). 
Indeed, given an arbitrary direction ft — (n x ,n y ,ri z ), 
parametrized in spherical coordinates (n x — sin x cos ■a/', 
n x = sin x sin ip, and n z = cosx), we can always find a ro- 
tation that brings |0) to \+)n and |1) to |— )^. Here \±)n 
are eigenvectors of the spin operator along the direction n, 
(na)|±)f2 = ±|±)^. Then, by performing a rotation before 
and after applying the measurement operators Eqs. (|1I2[) , 
we get the new measurement operators along ft, 

M&> = V^\+)nn(+\ + V^P\-)nn(-\, (10) 



M, 



(n) 



= Vq\+)nn{+\ + VPHnnH, 



(11) 



Similarly to Eqs. (01 H]) we get the effects E, 



(n) 



M^M™ and E$ 
bilities 



Mfo^ , as well as the proba- 



p(") 



MmWmWW), (12) 
(13) 



We can now calculate the Fisher information matrix, de- 
fined in general for a conditional probability distribution 
P(C|0) (where £ stands for the measurement outcomes 
and O a vectorial parameter) as the average Pq.q. = 
([a e ,lnP(C|e)] • [<9 0J lnP(C|e)]) over the values C (here 
don = d/dQi). The Fisher information is a measure of 
the sensitivity of our measurement scheme to the determi- 
nation of 9 and thus it is directly relevant for quantum 
tomography. In our case, we have O = (0,ip), £ € {m, fh}, 



giving a binary distribution pffi + Pf^' = !• This yields 



(H) _ T (n) _ 
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V 
(n) 



p (n) p (ff) 



[d e P^][d v P£\ (14) 
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0.0 



P 



rj(n) p (n) 
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3 (n) p(f?) 



d a P^ 



(15) 
(16) 



For the derivatives, we get 



doP^ = 



q-p 



sin 9 cos x — cos 9 sin x cos(ip — tp)] , 



d v Pm^ = q 2 P sin 9 sin x sin(<p - ijj) . 

These expressions, together with Eqs. (|T2| fT3]) . provide 
analytical expressions for calculating the Fisher matrix el- 
ements Eqs. (|14H16I) of a generalized measurement along 
any direction and for any state. As an example, in Fig. 
(JT|) we take x = 0, making all the elements of the Fisher 

(n) 

matrix zero except for J-g g , which is plotted as a function 
of p and q for 9 — tt/6. 

Generalized partial measurements are subject to the 
usual uncertainty principle: the precision in determin- 
ing one variable is limited by the Cramer-Rao bound 
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(a) 




Fig. 1: Plot of the function Tee for x 
function of p and q. 



and 6 = n/6 as a 



{e.g. (A9) 2 > 1/Tgg, if ip is known, or in general 



> T- 



where is the true value 



((e - e)(e - e) 4 ) 

of 6). 

In addition, another type of complementarity can be 
noticed. First, we see that no matter in which direction 
we perform the measurement, the four elements of the 
Fisher information matrix depend on (p — q) 2 (see e.g. 
Fig. [1]). Thus, the analysis based on the Fisher informa- 
tion confirms that no knowledge about the parameters 9, ip 
can be acquired near p = q. Therefore, to be able to do 
quantum tomography using generalized partial measure- 
ments it is preferable to have p and q as different from 
each other as possible (that is, one close to zero and the 
other one close to 1). But, on the other hand, we see that 
attempting to determine precisely the parameters charac- 
terizing the state increases the risk of not being able to 
reverse the measurement in either one of the two paths 
described above. Indeed, the likelihood of successful re- 
versal on a path to — > to is (1 — q)(l — p), while on the 
path to —> to is pq, and maximization of one leads to 
the other being zero. Now, since we have no knowledge 
about the state of the system, it is not possible to de- 
cide which path should be optimized; we can attempt to 
optimize them together by looking at the conditional en- 
tropy for successful measurements followed by reversals, 
-(l-p)(l-q)]n(l-p)(l-q)-pq]npq. At p = q = 1/2 
this quantity reaches its maximum value of In 2, and at 
this point the Fisher information is zero. 

Physical implementation. Physically, a POVM 
measurement can be regarded as a von Neumann mea- 
surement on another system (called "ancilla") which has 
previously interacted with the qubit, with the result de- 
scribed by a unitary U. The existence of this transforma- 
tion U is guaranteed by Naimark's dilation theorem but 
its realization in relation with experiment is usually not 
straightforward. Mathematically, given the basis |0,m), 



Fig. 2: (a) Implementation of the unitary U as a sequence 
of gates, (b) Schematic of a double-well implementation of 
the unitary U in a double-well potential. The vertical axis is 
an energy scale and the horizontal axis is the coordinate (a 
macroscopic superconducting phase difference, in the case of 
superconducting qubits). 



|0,to), |1, to), |1to), we can construct a unitary 
I+Z „ I-Z 



u 



® Up. 



(17) 



For example, for U q — ^1 — ql + i^fqY and U p = 
y/1 — pi + i^/pY we can check that the measurement op- 
erators M m and M™ are obtained by applying U and per- 
forming sharp (von Neumann) measurements with projec- 
tions | to) (to | and |to)(to|. Note also that even more gen- 
eral forms of the measurement operators can be obtained 
by using, instead of the UqS above, the most general form 
of a unitary on 2 x 2 matrices. 

To experimentally realize U, a straightforward imple- 
mentation is to simply consider two qubits and assume 
that one has available a universal set of gates. In Fig. 
H£a) we show the expansion of U in single-qubit X gates 
and two-qubit conditional gates. Such gates are currently 
available in many experimental setups, including super- 
conducting qubits (see e.g. [6] for the implementation of 
the conditional gates from Fig. [2](a) for flux qubits). 

Although possible in this way, it is still useful to look 
for an implementation via a single operation, which would 
save significantly on the time required to perform the mea- 
surement. The operation time can be an important re- 
striction, since many qubits have a limited coherence time. 
Such a physical implementation can be reached if we re- 
gard the ancilla as another degree of freedom of the same 
qubit. In this way, the number of qubits required is only 
one, which is significant if one weighs in the technical dif- 
ficulties of controlling two qubits and their coupling. Con- 
sider for example a symmetric double-well potential with 
adjustable barrier height (see Fig. [^b)). Initially, the 
barrier is high and the qubit is prepared in the left well in 
a superposition of the states |0) and |1) with energy differ- 
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ence hv, then it is lowered for some time r, allowing the 
qubit to tunnel, and then it is raised again. The tunnel 
matrix elements to and ti depend (via the WKB formula) 
on the energy (therefore on the state) of the qubit. We 
assume that the energy scales v 3> fo>ti are such that 
the lowering and raising of the barrier is adiabatic with 
respect to the qubit energy levels but instantaneous with 
respect to the tunneling rates. We refer to the left and 
right wells as m and m respectively. The Hamiltonian of 
the system during the time r is 

H = -^Z®I + t l-^^®X + t 1 l—^-®X. (18) 

We then solve the evolution problem in the frame rotating 
at the qubit frequency, i.e. we transform the wavefunction 
by —> exp(—ivZ ® I/2)t\iP), obtain the new Hamil- 
tonian, and evolve during r to finally get the evolution 
operator 

I + Z to 7 " to T 

U = — - — ® [cos -—I — i sin + 

I - Z . t\T T . . t\T 

H ® cos 1 — i sin X . 

2 L 2 2 J 

This obviously produces the required U with the identifi- 
cation q = sin 2 (to7"/2?i) and p = sin 2 (tir/2?i). 

Then, one should have available the projections |m)(m| 
and |m)(m|. In the field of superconducting qubits, only 
the first operator was available in standard switching- 
current measurements: the projection |m)(m| was real- 
ized if the measuring junction or SQUID did not switch 
in the running-wave state during the time r. If it did 
switch, the qubit was destroyed (similar to photon absorb- 
tion with optical qubits) due to quasiparticle generation 
at the point when the voltage reaches the superconducting 
gap. Recently, it was shown that QND measurements that 
could implement |m)(m| can be realized if the voltage is 
prevented to reach the gap by fast-feedback control elec- 
tronics [7]. Thus, with such devices it could be possible to 
implement the measurement operators described above. 

Conclusion. — We introduced a class of POVM 
measurements that generalizes the partial measure- 
ments demonstrated experimentally with superconducting 
qubits. We show that for these measurements it is possible 
to reverse (probabilistically) both of the measurement re- 
sults. We study the associated Fisher information and we 
propose a physical implementation of these measurements. 
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